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ABSTRACT 


The  generation  and  amplification  of  vortices  by  surface  inhomogeneities, 
both  in  the  form  of  surface  waviness  and  of  wall-normal  velocity,  is  investigated 
using  the  nonlinear  PSE  equations.  Transients  and  issues  of  algebraic  growth 
are  avoided  through  the  use  of  a  similarity  solution  as  initial  condition  for  the 
vortex. 

In  the  absence  of  curvature,  the  vortex  decays  as  y^l/x  when  flowing  over 
streamwise  aligned  riblets  of  constant  height,  and  grows  as  y/x  when  flowing 
over  a  corresponding  streamwise  aligned  variation  of  blowing/suction  transpi¬ 
ration  velocity.  However,  in  the  presence  of  wall  inhomogeneities  having  both 
streamwise  and  spanwise  periodicity,  the  growth  of  the  vortex  can  be  much 
harger.  In  the  presence  of  curvature,  the  vortex  develops  into  a  Gortler  vortex. 

The  “direct”  and  “indirect”  interaction  mechanisms  possible  in  wave-vortex 
interaction  are  presented.  The  “direct”  interaction  does  not  lead  to  strong 
resonance  with  the  flow  conditions  investigated.  The  “indirect”  interaction  leads 
to  K-type  transition. 
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1  Introduction 


We  present  results  from  an  investigation  into  the  laminar-turbulent  transition 
process  over  a  plate  with  wall  undulations,  wall  blowing  &  suction,  and  con¬ 
cave  curvature.  The  investigation  is  based  on  numerical  simulations  using  the 
nonlinear  PSE  equations.  The  curvature  is  sufficiently  small  to  allow  the  Bla- 
sius  profile  to  be  used  as  the  mean-flow  in  all  the  conditions  investigated.  The 
particular  topics  covered  in  this  report  are; 

1.  The  seeding  of  Gortler  vortices  by  streamwise  aligned  wall  corrugations 
and  blowing/suction  inhomogeneities. 

2.  The  receptivity  of  vortices  to  three  dimensional  wall  corrugations  and 
blowing/suction  inhomogeneities. 

3.  The  effect  of  curvature  on  three-dimensional  TS  waves. 

4.  Some  cases  of  nonlinear  wave-vortex  interactions. 

Two  geometries  are  considered.  Figure  la  displays  the  geometry  used  in 
the  study  of  topic  1.  The  wall  corrugations  vary  sinusoidally  in  the  spanwise 
direction,  z,  with  a  wavenumber  P,  and  have  the  crests  aligned  in  the  streamwise, 
or  X,  direction.  The  label  Xk  denotes  the  location  at  which  the  curvature  departs 
from  zero  and  increases  downstream  up  to  a  constant  final  value.  The  initial 
location  for  the  PSE  marching  solution,  xq,  is  upstream  of  Xfc. 

Figure  lb  displays  the  geometry  used  in  the  study  of  the  receptivity  of 
vortices  (Gortler  and  other  types)  to  wall  corrugations  which  vary  sinusoidally 
in  both  X  and  z  (topic  2).  This  receptivity  problem  is  steady  in  time.  Moreover, 
when  the  acoustic  source  is  added  (represented  in  the  figure  by  the  speaker),  the 
distributed  receptivity  mechanism  described  by  Crouch  [1]  is  activated,  creating 
Tollmien-Schlichting  waves.  These  waves  then  interact  with  the  vortex  (topic 

4). 

In  figure  1,  the  wall  undulation  can  be  interpreted  as  either  a  physical  un¬ 
dulation  present  on  the  plate  surface,  or,  in  the  case  of  blowing/suction,  as 
an  iso-level  of  wall-normal  velocity  on  a  smooth  plate.  The  height  of  the  wall 
undulation  is  of  the  order  of  0.1%  of  the  boundary  layer  thickness,  hence  the 
undulation  can  represent  imperfections  on  wing  surfaces  due  to,  for  example, 
temperature  or  direct  stresses  in  flight.  The  peak  surface  transpiration  veloc¬ 
ity  in  the  blowing  and  suction  study  is  of  order  of  1  x  10““*  t/oo.  in  agreement 
with  levels  used  in  Laminar  Flow  Control  experiments  (P.  Spalart,  personal 
communication) . 
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Figure  1:  Geometry  used  in  the  simulation.  Top,  the  corrugated  plate  with  an 
initial  flat  section,  followed  by  a  section  of  concave  curvature.  Bottom,  the  same 
geometry  as  above  with  wall  undulations  in  the  streamwise  direction  2ulded,  as 
well  as  acoustic  forcing. 
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The  main  results  corresponding  to  the  topics  listed  above  are: 

1 .  Wadl-aligned  corrugations  provide  an  efficient  means  for  the  generation  of 
Gortler  vortices. 

2.  Vortices  exhibit  growth  when  flowing  over  three  dimensional  wall  corru¬ 
gations  or  blowing/suction  inhomogeneities,  even  in  the  absence  of  cur¬ 
vature.  This  receptivity  persists  when  the  wavenumbers  of  the  wall  inho¬ 
mogeneities  differ  greatly  from  those  of  the  vortex,  i.e.  short  wavelength 
undulations  in  x  and  z  can  force  a  large  wavelength  vortex. 

3.  The  neutral  stability  curves  for  three-dimensional  TS  waves  enlarge  in  size 
as  curvature  is  increased,  and  eventually  extend  to  zero  frequency;  thus 
connecting  with  the  Gortler  vortices. 

4.  The  direct  and  indirect  wave-vortex  interaction  is  presented.  The  indirect 
interaction  leads  to  a  K-type  resonant  triad,  while  the  direct  interaction 
shows  no  wave-vortex  resonance  at  the  conditions  investigated. 

A  brief  description  of  the  PSE  equations  and  the  numerical  algorithm  em¬ 
ployed  to  solve  them  is  given  in  the  appendix  A.  We  also  describe,  in  appendix  B, 
the  construction  of  a  self-similar  solution  that,  when  used  as  initial  condition  for 
the  vortex  in  the  PSE  marching  procedure,  yields  results  free  of  transients  and 
of  algebraic  growth  (see  Schimd  and  Henningson  [2]).  The  self-similar  solution, 
thus,  allows  for  the  investigation  of  vortex  receptivity  to  distributed  forcing. 
Otherwise,  if  an  initial  condition  is  used  that  leads  to  strong  algebraic  growth, 
the  effects  of  the  distributed  forcing  will  be  masked.  A  similarity  solution  for 
the  Gortler  vortex  was  first  presented  by  Denier,  Hall  &  Seddougui  [3]  in  their 
investigation  of  vortex  forcing  from  a  plate  with  a  localized  hump. 

1.1  Reference  quantities,  and  geometry  description 

The  reference  length  chosen  is  6*  =  where  xJqoo  is  the  location 

at  which  the  Reynolds  number  =  VooX\qooIv*  equals  1  x  10®.  The  symbol 
*  denotes  a  dimensional  quantity.  This  reference  length  is  used  throughout 
this  report,  with  the  exception  of  appendix  B,  which  deals  with  the  self-similar 
solution. 

The  non-dimensional  curvature  is  /C  =  where  o^„,.  is  the  radius  of 

curvature.  Typical  values  of  AC  are  1, 2,4  and  8  x  lO""*.  In  the  case  of  a  subsonic 
wind-tunnel  test  with  Uoo  =  15  m/s,  the  reference  length  is  =  1.2  mm,  and 
the  radius  of  curvature  is  12,6,3  and  1.5  meters  respectively. 
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The  coordinate  along  the  streamwise  direction  is  x,  the  plate  normal  direc¬ 
tion  is  y,  and  the  spanwise  direction  is  z  The  corresponding  velocity  components 
are  u,  v,  and  w.  The  representation  of  the  flow-field  is  described  in  appendix 
A. 

All  modes  are  identified  with  the  three  integers  in  the  triplet  p  =  (/,n,  Jk), 
corresponding  to  the  indices  in  the  Fourier  expansion  6*"“®+**^*“*^*^*.  The  values 
of  a,  ^  and  u  are  prescribed  as  input  parameters  to  the  PSE  calculation.  The 
wavy-wall  surface  is  represented  by  the  function, 

f;  i;  +  C.C  (1) 

n=-N  k=-K 

where  “c.c”  denotes  complex  conjugate.  The  coefficients  W  are  complex  con¬ 
stants.  The  boundary  conditions  over  the  wavy  wall  are  transferred  to  y  =  0  via 
a  Taylor  series  expansion  (see  appendix  A).  The  blowing/suction  wall  normal 
velocity  is  described  by, 

E  E  5(0... +  c-c  (2) 

n=-Ar  k=-K 

The  boundary  conditions  at  the  wall  satisfied  by  the  disturbances  are  then. 

Up  =  0  Vp  =  S(o,n,k)  Wp  =  0  (3) 

The  homogeneous  boundary  condition  for  Up  and  Wp  implies  weak  suction 
rates  (Spalart  [4]).  Lastly,  in  the  presentation  of  the  results  we  use  the  local 
Reynolds  number  R  =  ^i/*x*/f/^  =  -v/fil,  and  the  nondimensional  frequency 
F  =  2irf*U^liy*,  where  /*  is  the  dimensional  frequency  in  Hz. 
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2  Validation  of  the  self-similar  vortex 

The  construction  of  the  self-similar  vortex  solution  is  described  in  appendix  B. 
Here,  we  present  the  results  which  prove  that  the  asumptions  and  approxima¬ 
tions  used  to  obtain  the  self-similar  formulation  are  sound. 

A  test  for  the  accuracy  of  the  self-similar  form  of  the  vortex  field  can  be 
made  by  running  the  PSE  using  as  initial  conditions  the  solution  to  (19)  and 
(20).  If  the  PSE  solution  remains  close  to  self-similar  over  a  downstream  range 
sufficiently  long  to  otherwise  display  algebraic  growth  and  decay,  then  the  self¬ 
similar  approximation  provides  an  acceptable  initial  condition  for  vortices  in 
both  direct  numerical  simulations  of  the  Navier-Stokes  equations,  and  in  PSE 
integrations.  The  test  was  run  using  the  linearized  PSE  equations  with  the 
initial  condition  given  at  /?  =  300.  The  riblets  on  the  wall  were  held  constant 
at  an  amplitude  of  0.004  while  in  the  suction  case,  the  level  wtis  held 

constant  at  1  x  10“®.  Figure  2  displays  the  maximum  amplitude  of  the  u  velocity 
component  as  function  of  Reynolds  number,  for  various  values  of  the  spanwise 
wavenumber  =  /3*^(jR=3oo)-  Note  the  abscence  of  transients  in  the  plot.  In 
the  case  of  the  wall  riblet,  the  amplitude  decreases  with  streamwise  distance 
due  to  the  following  reason:  the  maximum  of  the  u  velocity  occurs  at  j/  =  0 
where  the  boundary  condition  (21)  imposes  a  constant  amplitude  for  and 
the  physical  u  is  related  to  the  self-similar  by  the  factor  xq/x.  On  the 
other  hand,  constant  suction  leads  to  a  self-simialr  vortex  with  Umax  increasing 
proportional  to  y/x.  The  departure  of  the  PSE  solution  from  the  self-similar  one 
as  /?  is  increased  shows  the  error  introduced  by  freezing  the  coefficent  /^^(x/xo) 
in  the  self-similar  formulation. 

Figure  3  displays  the  velocities  u*e//(j?)  and  VseZ/C^)  at  the  initial  location 
R  =  300,  and  downstream,  at  R  =  1136  and  1670.  The  dots  trace  the  shape  of 
the  self-similar  functions  u,eZ/(>/)  and  v,e//(?/),  while  the  solid  and  dashed  lines 
show  the  corresponding  profiles  obtained  by  extracting  the  u,e//  and  profiles 
from  those  given  by  the  PSE  equations  at  R  =  300  to  R  =  980  and  R  =  1600. 
Both  the  u  and  the  v  velocity  components  computed  with  the  PSE  collapse  close 
to  the  self-similar  profiles.  This  agreement,  along  with  the  abscence  of  algebraic 
growth  or  other  transients,  validates  the  assumptions  and  approximations  used 
in  generating  the  self-similar  solution. 
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Figure  3:  Self-Similar  profiles 


3  Vortex  Seeding 

The  generation  of  Gortler  vortices  by  a  localized  hump  was  investigated  by 
Denier,  Hall  &  Seddougui  [3],  and  found  that  the  mechanism  was  inefficient 
in  generating  vortices.  Later,  Bassom  and  Hall  [5]  corrected  some  errors  and 
extended  the  work,  and  found  that  distributed  forcing  functions  are  much  more 
efficient  in  the  generation  of  Gortler  vortices. 

We  study  the  birth  and  growth  of  the  Gortler  vortices  that  develop  on  top 
of  the  corrugated  (i.e.  ribbed)  plate  with  a  concave  curvature.  Initially,  the 
curvature  of  the  corrugated  plate  is  zero,  and  is  increased  smoothly  from  zero 
at  /?  =  300  to  the  maximum  constant  value  &t  R  >  500  with  a  half-period 
cosine  function. 

At  the  onset  of  curvature,  the  flow  is  composed  of  the  Blasius  flow  plus  the 
self  similar  solution  discussed  above.  Downstream,  Gortler  vortices  develop  with 
the  same  spanwise  wavenumber  as  the  wall  corrugation.  The  self-similar  solution 
has  a  non-zero  projection  onto  the  Gortler  eigenmode,  and,  thus,  provides  the 
initial  amplitude  for  the  vortex. 

The  five  c<ises  studied  are  categorized  by  the  variation  in  x  of  the  wall 
corrugations.  The  case  ‘infinite’  has  constant  amplitude  for  all  R  >  300;  the 
case  ‘400-500’  has  a  constant  amplitude  up  to  /?  =  400,  then  decreases  smoothly 
to  zero  between  400  <  R  <  500  and  remains  zero  afterwards;  the  cases  ‘400-600’ 
and  ‘400-800’  and  ‘300-400’  are  of  equal  nomenclature.  For  each  case,  5  runs  are 
made  with  the  curvature  taking  on  values  of  /C  =  0, 1  x  10“^, 2  x  10“^, 4  x  10"“* 
and  8x  10“^.  The  spanwise  wavenumber  in  all  cases  is  =  0.15.  The  calculation 
is  maintained  purely  linear  by  artificially  neglecting  the  nonlinear  terms  that 
arise  from  a  finite  amplitude  vortex. 

Figure  4  displays  the  maximum  amplitude  of  the  u  component  of  velocity  as 
function  of  Reynolds  number.  The  plots  on  the  left  column  of  figure  1  display 
the  Umax  amplitude  of  the  vortex  mode  versus  Reynolds  number  for  the  cases 
‘infinite’,  ‘400-600’,  and  ‘300-400’.  The  Gortler  vortex  amplitude  is  ploted  on 
the  right  column.  This  amplitude  is  computed  by  subtracting  from  the  flow-field 
the  flow  field  due  to  the  wall  corrugation  in  the  absence  of  curvature. 

The  ‘infinite’  case  produces  the  fastest  growth  of  the  Gortler  mode.  This 
fact  suggests  that  the  vortex  is  receptive  to  the  forcing  from  the  wall  corrugation 
over  an  extended  streamwise  extent.  Indeed,  the  extend  seems  larger  than  the 
region  of  strong  nonlocal  receptivity  (Crouch[l],  Crouch  &  Bertolotti  [6])  for 
traveling  waves.  The  growth  rates  in  all  cases  asymptote  to  a  single  curvature- 
dependent  value  aX  R  =  1400:  the  difference  in  amplitude  at  this  position  from 
case  to  case  is  due  to  the  different  corrugation  geometry. 


8 


The  ‘kink’  in  the  amplitude  curves  of  the  ‘400-600’  and  ‘300-400’  cases  occurs 
at  the  streamwise  location  where  the  decaying  wall-mode’s  Umar  amplitude, 
which  is  maximum  at  the  wall,  equals  the  growing  Gortler  amplitude,  which 
has  a  maximum  at  about  y  =  2.5. 

An  interesting  phenomena  can  be  seen  in  the  amplitude  plots  of  the  Gortler 
vortex  in  the  ‘300-400’  case.  The  kink  in  the  amplitude  curve  in  the  neighbor¬ 
hood  oi  R  =  800  divides  the  upstream  region  dominated  by  the  wall  mode  from 
the  downstream  region  dominated  by  the  Gortler  vortex.  The  streamwise  po¬ 
sition  of  this  kink  shifts  downstream  when  curvature  is  increased.  This  suggest 
that  the  wall-mode  retards  the  growth  of  the  Gortler  vortex. 

4  Vortex  Receptivity 

We  study  the  effect  of  surface  roughness  on  the  growth  of  vortices.  I  have 
used  the  word  “seeding”  in  the  previous  section  and  the  word  “receptivity”  in 
this  section  to  differentiate  between  the  types  of  forcing  given  to  the  vortex. 
Seeding  occurs  when  the  forcing  comes  directly  from  the  boundary  condition, 
and  hence,  is  a  linear  process  (the  word  “seeding”  was  suggested  by  M.  Morkovin 
in  a  private  communication).  Receptivity  occurs  when  the  forcing  comes  from 
the  nonlinear  interaction  between  other  modes.  In  the  absence  of  nonlinearity 
(or  in  the  limit  of  small  amplitude),  these  modes  do  not  affect  the  vortex  at 
all.  In  our  case  these  modes  are  due  to  wall  undulations  or  wall  suction  having 
wavenumbers  different  from  that  of  the  vortex. 

In  order  to  isolate  the  receptivity  effect  from  centrifugal  instabilities,  we 
consider  the  case  of  a  flat  plate  at  zero  angle  of  attack.  The  plate’s  surface 
contains  a  spanwise  periodic  rib  with  streawise  aligned  peaks  and  valleys.  This 
rib  generates  the  vortex  -  and  at  the  initial  marching  location  the  vortex  form 
is  given  by  the  self-similar  solution.  The  height  of  the  rib  is  constant  in  x,  and 
corresponds  to  the  coefficient  W(o,o,2)  =  1  x  10~^  (see  equation  1).  In  addition 
to  this  rib,  we  add  the  source  for  other  steady  disturbances,  either  in  the  way  of 
additional  surface  undulations,  or  using  blowing  and  suction  periodic  in  x  and 
z.  The  height  or  strength  of  this  additional  source  is  increased  from  0  to  its  final 
value  using  a  ramp  function  starting  aX  R  =  425  and  ending  at  A  =  700.  The 
gradual  introduction  of  these  modes  reduces  the  amount  of  algebraic  growth 
that  is  introduced  into  the  evolution  of  the  vortex. 

The  first  set  of  results  are  for  the  case  in  which  the  wall  contains  one  addi¬ 
tional  Fourier  mode  mode  (besides  the  rib),  having  a  streamwise  wavenumber  a 
that  is  varied  parametrically,  and  a  spanwise  wavenumber  0  that  is  half  that  of 
the  vortex.  The  indices  for  the  wall  mode  are  of  the  form  (0,n,  1),  while  that  of 


the  vortex  is  (0,0, '2).  Thus,  the  interaction  of  mode  (0,n,  1)  with  its  companion 
mode  (0,  — n,  1)  (i.e.  from  symmetry  in  0)  creates  a  forcing  at  (0,0,2). 

Figure  5  displays  the  results  for  the  case  of  a  wavy  wall  with  Fourier  coef¬ 
ficient  W(o,n,i)  =  2  X  10“^,  and  the  case  of  blowing/suction  with  a  coefficient 
5(0, n.i)  =  5  X  10“®,  for  different  values  of  a.  In  the  wavy- wall  case,  increasing 
a  increases  the  forcing  efficiency.  For  a  —  1.0  the  vortex  increases  one  order  in 
magnitude.  In  the  suction  case,  the  opposite  trend  exists;  reducing  a  increases 
^he  forcing.  The  explanation  for  this  opposite  trends  can  be  seen  in  the  velocity 
profiles  of  the  steady  modes,  shown  in  figure  6.  Shortening  the  wavelength  of 
the  wavy- wall  mode  increases  the  slope  dufdy  at  the  wall,  which  is  multiplied 
by  the  wall  undulation  height  to  produce  the  forcing  (see  equation  (11)).  On 
the  other  hand  decreasing  a,  while  holding  the  suction  level  constant,  increases 
the  u  component  of  the  velocity,  hence  the  forcing  in  the  interior  of  the  domain. 

The  second  set  of  results  we  present  contain  a  pair  of  steady  modes  having 
a  difference  in  spanwise  wavenumber  equal  to  the  vortex’s  wavenumber, 

W(0,n,m)  =  2  X  10'^  W(0,n.m-(.2)  =  *  2  X  10’^ 

<5(0, n,m)  =  5  X  10  ^  <5(0,n,m+2)  =  t  5  X  10  ^ 

The  quadratic  nonlinearity  then  produces  a  forcing  at  (0,0,2)  for  any  value 
of  n  and  m.  In  this  way,  short  wavelength  disturbances  can  directly  influence 
long  wave-length  disturbances.  This  type  of  coupling  was  used  by  Crouch  [7] 
in  a  study  of  TS  wave  receptivity  to  short-scale  waviness.  The  pheise-shift  of 
90  degrees  between  the  two  steady  modes  leads  to  a  more  effective  forcing  than 
the  case  were  both  Fourier  coefficients  W  are  purely  real. 

Figure  7  displays  the  result  for  the  combinations  n  =  l,m  =  2,  and  n  = 
l,m  =  6,  with  a  =  1.0  and  ^  =  0.15.  The  vortex  has  a  spanwise  wavenumber 
oi  0  =  0.30,  as  in  the  previous  case.  For  reference,  the  case  0  +  0  form  figure  5 
is  included  in  the  plot.  The  large  vortex  amplitude  created  by  the  “difference 
interactions”  shows  that  the  coupling  between  short  and  long  wavelengths  can 
be  strong. 

Figure  8a  show  the  effect  of  the  surface’s  undulation  height  on  the  evolution 
of  the  vortex.  A  height  as  small  as  2.5  x  10“^  can  produce  a  noticeable  increase 
in  vortex  amplitude.  For  reference,  it  is  useful  to  consider  dimensional  quan¬ 
tities.  In  a  slow-speed  wind-tunnel  with  a  free-stream  velocity  of  10  m/s,  the 
reference  length  is  So  =  1.5  mm  (i.e.  at  /?  =  1000),  thus  the  surface  undula¬ 
tions  are  roughly  15  micrometers.  This  small  value  suggest  that  the  distributed 
forcing  form  small  plate  undulations  can  introduce  vortices  in  cases  where  the 
curvature  is  nominally  zero.  In  particular,  distributed  forcing  can  be  a  cause 
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Wavy-Wall 
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Figure  5:  Umax  amplitude  versus  Reynolds  number  of  a  vortex  with  spanwise 
wavenumber  2^  flowing  over  a  surface  with  wavenumber  Top,  surface 

mode  corresponds  to  surface  undulation;  bottom,  surface  mode  corresponds  to 
blowing/suction.  The  line  labeled  “rib  only”  shows  the  amplitude  of  the  vortex 
flowing  over  the  streamwise  ribs  only. 
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Figure  6:  u  velocity  profiles  at  /Z  =  600  due  to  a  wavy- wall,  left,  and  suction, 
right,  with  spanwise  wavenumber  ^  ~  0.3  and  different  values  of  streamwise 
wavenumber  a. 
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Figure  7;  Umax  amplitude  versus  Reynolds  number  of  a  vortex  with  spanwise 
wavenumber  2/?  flowing  over  a  surface  with  modes  (0,a,n/?),  where  a  =  1.0, 
and  =  0.15.  Top,  surface  modes  corresponds  to  surface  undulation;  bottom, 
surface  mode  corresponds  to  blowing/suction.  The  line  labeled  “rib  only”  shows 
the  amplitude  of  the  vortex  flowing  over  the  streamwise  ribs  only. 
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of  the  Klebanoff  mode  [8],  in  addition  to  the  possibility  of  algebraic  growth  as 
investigated  by  Herbert  and  Lin  [9]. 

Figure  8b  shows  the  velocity  profiles  of  the  vortex  for  the  case  n  =  l,m  =  1 
of  figure  8a.  The  u  profiles  agrees  in  shape  with  the  experimental  measurements 
of  Kendall  [10]  of  the  Klebanoff  mode. 
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Figure  8:  (a)  EflFect  of  surface  height,  |W|,  on  the  Umax  amplitude  for  a  vor¬ 
tex  with  spanwise  wavenumber  flowing  over  a  surface  with  wavenumbers 
(0,  Of,  4^)  and  (0,  a,  2^).  (b)  vortex  «  and  v  velocity  profiles  at  various  stream- 
wise  locations  for  the  case  |W|  =  5  x  10“^. 


5  The  effect  of  curvature  on  three-dimensional 
TS  waves 

We  turn  our  attention  now  to  three-dimensional  TS  waves  and  the  effect  that 
curvature  has  on  them.  The  motivation  for  this  study  is  two-fold.  First,  we  look 
at  the  effect  of  curvature  on  disturbances  that  have  a  wave-vector  orientation 
in  between  that  of  vortices  (i.e.  purely  along  2)  and  that  of  two-dimensional  TS 
waves  (i.e.  purely  along  x).  Second,  the  waves  studied  will  be  candidates  for 
the  nonlinear  wave-vortex  study,  discussed  in  the  next  section. 

In  this  study,  the  wall  corrugations  are  omitted,  leaving  a  smooth  plate. 
The  boundary  conditions  are  homogeneous,  leading  to  the  classical  eigenvalue 
problem  for  the  TS  waves.  We  solve  for  an  augmented  eigenvector  composed 
of  {u,v,du/dx,dv/dx}  and  an  augmented  eigenvalue  (a,dafdx)  using  the  “lo¬ 
cal  procedure”  described  in  Bertolotti,  Herbert  &  Spalart  [11].  This  procedure 
yields  an  improved  solution  over  the  standard  Orr-Sommerfeld  solution:  in  par¬ 
ticular,  it  captures  the  increased  effect  of  non-parallelism  on  three-dimensional 
waves  (Bertolotti  [12]).  However,  the  procedure  becomes  invalid  in  the  limit  of 
low-frequency  and  long-wavelengths,  where  the  governing  equation  cannot  be 
reduced  to  ordinary  differential  form.  In  particular,  the  validity  of  the  local  pro¬ 
cedure  breaks  down  when  the  streamwise  wavelength  of  the  eigenmode  becomes 
of  order  0(R),  hence  the  procedure  cannot  reach  the  limit  of  pure  vortices. 

We  use  the  local  procedure  because  it  allows  the  neutral  curves  to  be  traced 
efficiently  in  the  R-F  plane.  Using  the  PSE  solution,  on  the  other  hand,  requires 
many  runs  at  different  frequencies  in  order  to  “raster-scan”  the  neutral  curve, 
like,  if  you  will,  the  electron  beam  on  a  TV  screen. 

The  neutral  stability  curves,  based  on  udy,  for  TS  waves  having  a  fixed 
spanwise  wavenumber,  /?  =  0.30  and  =  1.0  are  shown  in  figure  9.  In  the 
absence  of  curvature  the  low-frequency  TS  waves  are  stable,  but  with  some 
curvature  present  a  second  region  of  instability  appears  near  F  =  0.  Focusing 
on  the  =  0.3  case,  we  see  that  as  the  curvature  increases,  two  neutral  curves 
form:  the  upper  is  the  continuation  of  the  curve  for  the  flat  plate;  the  second, 
at  low  frequencies,  displays  the  fact  that  centrifugal  instability  is  amplifying  the 
Gortler  vortices  (F  =  0)  as  well  as  modes  with  low  frequency.  At  a  curvature 
value  of  8  X 10““*  the  two  curves  merge.  At  this  point  there  is  a  smooth  connection 
between  TS  waves  and  the  Gortier  vortices. 

The  numbers  along  the  outer  neutral  curve  represent  the  wave-angle  of  the 
TS  wave  when  AC  =  8  x  10"*.  The  table  below  gives  the  phase-speed  of  the  TS 
wave  at  the  location  where  the  symbol  of  the  wave  angle  appears  in  figure  9. 
Note  the  increase  in  phase-speed  as  the  frequency  is  lowered,  and  the  wave-angle 
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Figure  9:  Effect  of  curvature,  k,  on  the  neutral  stability  curves  for  3-D  Tollmien- 
Schlichting  waves  over  a  smooth  plate.  The  numbers  along  the  «  =  8  x  10"^ 
curve  show  the  angle  (deg.)  of  the  wave-vector  to  the  free-stream  direction. 
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is  increased. 


Neutral  curve,  ^  ~  0.3,  k  =  8  x  10  * 


Wave  Angle  (deg.) 

Phase  Speed 

34.4° 

0.4016 

43.2° 

0.3984 

44.8° 

0.3930 

53.8° 

0.4019 

59.1° 

0.4014 

65.4° 

0.4174 

69.4° 

0.4226 

72.8° 

0.4401 

83.5° 

0.5114 

85.4° 

0.5735 

Neutral  curve,  = 

1.0,  K  =  8  X  lO-'* 

Wave  Angle  (deg.) 

Phase  Speed 

85.1 

0.6180 

86.2 

0.6510 

86.5 

0.6424 

86.7 

0.6829 

TS  waves  with  a  ^  =  1.0  have  a  1/3  shorter  spanwise  wavelength  than  those 
with  P  =  0.3,  and,  hence,  are  more  densely  packed.  In  the  Blasius  boundary 
layer,  these  waves  are  stable.  Furthermore,  the  well-known  types  of  nonlinear 
wave  interactions,  i.e.  K-type  and  H-type,  involve  waves  with  /?  roughly  equal 
to  a,  hence  in  the  range  j3  <  0.5.  Thus,  we  must  look  at  the  growth  rates  to 
make  some  comments  on  importance  of  these  densely  packed  waves. 

Figure  10  compares  the  growth  rates  of  a  vortex  with  /S  =  0.3  and  those  of 
a  TS  wave  with  ^  =  0.3  and  F  =  20.  The  vortex  grows  over  a  fully  corrugated 
wall,  while  the  waves  grow  over  a  smooth  wall.  The  growth  rate  of  the  vortex 
is  greater  than  that  of  the  TS  wave,  being  roughly  twice  as  large  for  the  case 
of  curvature  =  k  =  8  x  10“^.  For  comparison,  the  maximum  growth  rate  of  a 
2-D  TS  wave  with  F  =  60  is  7  =  0.0141.  Thus,  the  low-frequency  3-D  TS  wave 
undergo  slower  growth,  but  do  so  over  a  more  extended  streamwise  distance. 

We  conclude  this  section  with  figure  11,  which  shows  the  variation  of  the 
receptivity  coefficient  in  Crouch’s  receptivity  model  [1]  with  curvature  for  a  2-D 
TS  wave  at  F  =  60  and  a  3-D  TS  wave  with  F  =  60  and  ^  =  0.15.  The 
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Figure  10:  Comparison  of  growth  rates,  based  on  Umax,  of  a  Gortler  vortex  with 
0  =  0.3,  and  a  TS  wave  of  frequency  F  =  20  and  ^  =  0.3,  for  various  values  of 
curvature. 
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Figure  1 1 :  Effect  of  curvature  on  the  receptivity  coefficient  for  a  two-dimensional 
TS  wave,  /9  =  0,  and  a  three-dimensional  TS  wave,  ^  =  0.15  with  frequency 
F  =  60. 

data  shows  a  weak  influence  of  curvature  on  the  receptivity  mechanism  (at  this 
frequency). 
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6 


Some  case-studies  of  wave-vortex  interac¬ 
tion. 

The  interaction  between  vortices  and  traveling  waves  in  boundary-layers  can 
play  a  major  role  in  determining  the  transition  location  on  swept  wings.  Partly 
for  this  reason,  this  interaction  has  been  the  focus  of  numerous  studies. 


The  two  possible  ways  a  vortex  and  a  wave  can  inter2u:t  are  shown  in  the 
figure  above.  The  left  diagram  shows  the  triad  interaction  between  a  vortex 
and  two  oblique  TS  waves  of  equal  frequency  and  streamwise  wavenumber, 
but  opposite  spanwise  wavenumber.  We  shall  refer  to  this  interaction  as  the 
“direct  interaction”,  since  in  a  perturbation  approach  employing  an  expansion 
in  amplitude,  the  waves  interact  at  first  order.  The  right  diagram  shows  a  triad 
interaction  between  a  vortex,  a  two-dimensional  traveling  wave  and  a  three- 
dimensional  traveling  wave.  The  traveling  waves  have  equal  frequency  and  a 
nearly  equal  streamwise  wavenumber.  The  oblique  vector  is  of  a  lighter  shade  of 
gray  to  indicate  that  its  presence  is  not  necessary  at  the  onset  of  the  interaction 
-  this  wave  will  be  generated  by  the  interaction  between  the  vortex  and  the  2D 
wave,  and  will  rise  in  amplitude  to  close  the  triad  inter2w:tion.  We  shall  refer  to 
this  interaction  as  the  “indirect  interaction”,  since  in  a  perturbation  approach 
employing  an  expansion  in  amplitude,  the  waves  interact,  initially,  at  second 
order. 

The  direct  interaction  has  been  investigated  more  thoroughly  than  the  indi¬ 
rect  interaction.  Studies  of  the  direct  interaction  include  those  of  Hall  and  Smith 
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[13]  who  employed  a  triple-deck  asymptotic  expansion  in  an  analysis  entitled 
“on  strongly  nonlinear  vortex/wave  interactions  in  boundary- layer  transition”; 
Davis  and  Smith  [14],  who  extended  the  investigation  to  three-dimensional 
boundary  layers;  Thumm,  Wolz  and  Fasel  [15],  who  looked  at  the  interaction  in 
a  compressible  boundary-layer  using  a  direct  numerical  simulation,  starting  with 
only  the  oblique  waves;  Chang  and  Malik  [16],  who  looked  at  the  same  case  using 
the  PSE  equations;  and  Spalart  [4],  who  employed  a  direct  numerical  simulation 
to  look  at  the  generation  of  cross-flow  vortices  by  suction  non-uniformities,  and 
the  subsequent  breakdown  induced  by  interaction  with  a  TS  wave. 

With  the  exception  of  Spalart  [4],  the  vortex  in  these  investigations  was 
generated  by  the  nonlinear  interaction  of  two  oblique  TS  waves.  Through  this 
process,  the  vortex’s  motion  is  “in  tune”  with  the  forcing  provided  by  the  trav¬ 
eling  waves,  and  the  triad  interaction  is  optimal  from  a  phase  point  of  view.  In 
our  investigation  the  vortex  and  the  traveling  waves  exist  independently  at  the 
initial  marching  location,  and  may  not  be  as  much  “in  tune”.  This  difference 
may  explain  the  lack  of  resonance  found  in  the  results  presented  below. 

The  indirect  interaction  was  studied  by  Bertolotti  [12].  The  PSE  equations 
where  employed  in  a  parametric  study  of  the  effect  of  initial  amplitude  levels 
and  vortex  spanwise  wave-number  on  the  location  of  “transition”,  (defined  as 
the  onset  of  rapid  spectrum  filling).  It  v/as  shown  that  vortex  and  TS  wave 
behaved  linearly  up  to  the  streamwise  location  where  the  the  oblique  wave 
matched  amplitudes  with  the  vortex.  Downstream  of  this  location  a  strong 
resonance  ensued,  wherein  the  growth  rates  of  the  vortex  and  of  the  oblique 
wave  where  equal  to  one  another,  and  were  one  order  of  magnitude  larger  than 
those  given  by  linear  calculations.  These  growth  rates  matched  those  given  by 
Floquet  theory  for  K-type  transition.  Figure  12,  taken  from  Bertolotti  [12], 
displays  these  results.  The  initial  conditions  consists  of  the  (2,0)  TS  mode,  and 
the  (0,1)  vortex.  (In  this  work  the  disturbances  were  phased-locked,  thus  the 
index  (2,0)  stands  for  (2,2,0)  in  the  current  terminology).  The  dots  represent 
the  amplitude  of  the  K-type  secondary-instability  wave  with  equal  parameters 

F,  and  Ats  as  in  the  vortex  case.  The  agreement  in  slope  between  the 
dots  and  growth- rates  of  the  vortex  and  the  (2,1)  mode  inuicates  that  these 
two  modes  undergo  resonance  as  in  the  Floquet  model.  This  figure  was  not 
published  because  the  investigators  felt  that  the  way  in  which  the  vortex  was 
initialized,  namely  by  taking  the  limit  of  F  — >  0  holding  0  constant  in  the  local 
procedure  was  too  arbitrary  a  choice. 

In  our  study  of  wave- vortex  interaction,  we  will  include  the  receptivity  mech¬ 
anism.  In  this  way,  the  initial  conditions  are  specified  by  a  small  set  of  data, 
namely  the  description  of  the  wall  geometry  and  of  the  acoustic  field  present 
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Figure  12:  Wave-vortex  indirect  interaction:  Amplitudes  of  selected  modes  ver¬ 
sus  Reynolds  number.  Initial  conditions  specified  for  the  vortex  and  the  (2,0) 
TS  wave  at  frequency  F  =  30.  The  dots  show  the  results  from  K-type  secondary 
instability  given  by  Floquet  theory.  (Plot  taken  from  Bertolotti’s  PhD  Thesis, 
page  165) 
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in  the  free-stream.  We  employ  Crouch’s  model  [1]  for  the  generation  of  the  TS 
waves  form  the  interaction  of  the  acoustic  field  and  the  wall  geometry.  (See 
Bertolotti  and  Crouch  [17]  for  details  on  the  PSE  implementation.) 

Two  TS  frequencies  were  investigated,  one  at  F  =  60,  representing  the  higher 
frequency  modes,  the  other  at  F  =  20,  representing  the  lower  frequencies.  As 
shown  in  figure  9,  low-frequency  TS  waves  are  amplified  when  concave  curvature 
is  present,  and,  thus,  might  play  an  important  role  in  the  transition  process 
over  concave  surfaces.  The  curvature  is  held  fixed  at  k  =  2  x  10"^.  This  small 
value  of  curvature  is  chosen  because  it  leads  to  Gortler  vortices  with  growth 
rates  similar  to  those  of  the  TS  waves.  If  the  curvature  is  increased  to  higher 
values,  say  8  x  10“'^,  then  the  growth  of  the  vortex  is  too  rapid  to  allow  for 
significant  interaction  with  the  traveling  waves,  and  at  large  amplitudes  the 
vortex  generates  the  transitional  state  suggested  and  investigated  by  Hall  and 
Horseman  [18],  wherein  the  inflectional  mean-flow  profile  induced  by  the  vortex 
gives  rise  to  inviscid  (i.e.  Rayleigh)  instabilities. 

At  each  frequency,  two  selections  of  surface  roughness  are  used,  a  lower  one 
with  peak-to-peak  variations  in  the  range  of  1  x  10~^^o»  and  a  higher  one  with 
variations  in  the  range  of  1  x  10“^^o. 

The  parameters  describing  each  case  are  listed  in  the  table  below: 

Direct  Interaction 


Indirect  Interaction 


Since  the  mean-flow  is  independent  of  2,  symmetric  disturbances  in  z  are 
assumed.  To  obtain  peak-to-peak  values,  one  must  multiply  W(o,n,o)  and  W(o,o,n) 
by  2,  and  W(o,n,m)  by  4  to  take  into  account  the  complex  conjugate  modes. 
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The  amplitude  evolution  of  the  vortex,  the  traveling  waves,  and  the  mean- 
flow  distortion  are  displayed  in  figure  13  for  the  frequency  F  =  20  and  in  figure 
14  for  the  frequency  F  =  60.  (These  frequencies  values  refer  to  the  traveling 
waves).  Both  the  low  and  high  wall  roughness  cases  lead  to  a  vortex- wave 
resonance  in  the  case  of  the  indirect  interaction.  The  location  of  this  resonance 
is  only  slightly  affected  by  the  level  of  the  wall  undulations. 

On  the  other  hand,  the  cases  involving  the  direct  interaction  do  not  lead  to 
a  wave-vortex  interaction.  The  traveling  wave  (i.e.  mode  (2,2,1))  is  dampened 
past  R  =  1200,  most  likely  due  to  the  presence  of  the  large  mean-flow  distortion 
induced  by  the  vortex.  The  higher  growth  rate  of  the  vortex  in  the  “High”  case 
is  due  to  the  receptivity  of  the  vortex  to  the  (2,2,1)  wall  mode,  as  discused 
in  section  5.  The  absence  of  a  wave-vortex  interaction  is  surprising,  since  the 
papers  cited  at  the  beginning  of  this  section  predict  a  strong  resonance.  A 
possible  reason  was  thought  to  be  due  to  an  unfavorable  phase  relation  between 
the  vortex  and  the  3-D  TS  wave.  However,  changing  the  phase  by  90  and  180 
degrees  did  not  change  qualitatively  the  results.  Thus,  the  reason  for  the  lack 
of  resonance  remains  at  the  moment  unresolved. 

The  surface  heights  considered  here  are  of  the  same  order  of  magnitude  as 
adhesive  tape.  The  indirect  interaction  calculations,  then,  can  help  understand 
in  a  quantitative  way  Klebanov’s  experiment  [19],  in  which  strips  of  adhesive 
tape  where  placed  under  the  vibrating  ribbon  to  help  steady  the  transitional 
flow  structure. 


7  Conclusions 

The  response  of  vortices  to  surface  inhomogeneities,  both  in  the  form  of  surface 
waviness  and  of  wall-normal  velocity,  has  been  successfully  investigated  using 
the  nonlinear  PSE  equations.  Transients,  and  issues  of  algebraic  growth,  have 
been  avoided  through  the  use  of  a  similarity  solution  as  initial  condition  for  the 
vortex. 

In  the  absence  of  curvature,  the  vortex  decays  as  yj\jx  when  flowing  over 
streamwise  aligned  riblets  of  constant  height,  and  grows  as  y/x  when  flowing 
over  a  corresponding  streamwise  aligned  variation  of  blowing/suction  at  the 
wall.  However,  in  the  presence  of  wall  inhomogeneities  having  both  streamwise 
and  spanwise  periodicity,  the  growth  of  the  vortex  can  be  much  larger.  In  the 
presence  of  curvature,  the  vortex  develops  into  a  Gortler  vortex. 

The  “direct”  and  “indirect”  interaction  mechanisms  possible  in  wave- vortex 
interaction  are  presented.  The  “direct”  interaction  does  not  lead  to  strong 
resonance  in  the  flow  conditions  investigated.  The  “indirect”  interaction  leads 
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Figure  14:  «mox  amplitudes  versus  Reynolds  number  with  a  traveling-wave  fre¬ 
quency  of  F  =  60. 
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to  K-type  transition. 
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9  Appendix  A:  The  governing  equations 

We  eliminate  pressure  by  taking  the  curl  of  the  Navier-Stokes  equation.  Since 
the  flow  field  is  harmonic  in  2,  we  can  eliminate  w  using  continuity,  and  then 
take  the  combinations  of  d/dx{k  component)  -  d/dz{\  component)  ,  and  5/9z(j 
component)  of  the  vorticity  equation  to  obtain  two  equations  governing  u  and 

V. 

The  resulting  formulation  offers  two  advantages  over  the  original  equations 
for  the  full  set  of  primitive  variables;  reduced  computational  work,  and  PSE 
equations  that  are  free  from  numerical  instabilities  at  short  marching  steps. 

We  separate  the  flow  filed  vt  into  the  basic-flow,  Vb(i,j/)  and  the  remain¬ 
der  v(x,  j/,  z,  f),  which  we  call  the  disturbance.  Our  basic-flow  is  the  Blasius 
boundary  layer,  and  terms  involving  the  second  derivative  with  respect  to  x  are 
of  negligible  magnitude.  Excluding  these  terms  the  governing  equations  for  the 
disturbance  components  u  and  v  are: 

,  JT  ^  ^  dZs  ,^du  dv , 


dZadv  d^Zs  dUB,dH 
dy  dx  ^  dxdy^  ^  dx 


d^v  d^u 

dz'^  dxdy 


f  d^Vg  dVB  d^u  d^u  \ 

{dxdy^''^  dx^dz^  dy^’j 

4.2ic\u  +  + 

®  dx"^  dz^  dxdy  dx  dy 

dUgdv  dVB.du  d^u  _ 

dy  dx  dx  dy  dx  ^  dxdy 

,.d  d  ^d.  ,  d  d  d  .,dC  d(. 

-% + ’'s; + 

^  9.  , 

dz  dx  dz  dy  dz  dz  dx  dx  dx  dy  dx  dz  ^ 

du  d  dv  d  dw  d  du  d  dv  d  dwd 

dz  dx^  dz  dy  dz  dz  'dxdx  dx  dy  dx  dz 


,d  ..  d  d  1 


d^u  d^v  . 
dz"^  dxdy 
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(46) 


dy"^  dz  dz^  dz  dx  dz  dy'^  dz  dz 
d  d  d  .dr}  du  d  dv  d  Ow  d  . 

^  ^dy  ^  ^ dz^ dz  ^dzdx  ^  dz  dy  ^  dz  dz  ^ 

We  have  placed  the  nonlinear  terms  involving  the  product  of  disturbance  quan¬ 
tities  as  right-hand-side  terms  because  these  terms  will  be  lagged  in  the  iterative 
process  used  to  solve  the  system  of  equations. 

The  basic-flow  terms  on  the  third  line  of  equation  (4a)  are  negligible  within 
boundary-laye  theory.  When  the  disturbance  has  velocities  components  of  equal 
magnitude,  such  as  a  traveling  wave,  then  these  terms  are  negligible.  However, 
when  the  disturbance  has  u  and  v  velocity  components  that  scale  like  those  of 
the  Blasius  flow,  such  as  vortices,  then  the  terms  on  the  third  line  of  (4a)  become 
of  equal  order  to  terms  in  the  first  two.  Thus,  in  cases  where  the  disturbance 
field  includes  both  waves  and  vortices,  as  in  the  present  study,  these  terms 
should  be  kept.  In  passing  we  note  that,  for  steady  vortices,  equations  (4a)  and 
(4b)  equal  Hall’s  equations  for  Gortler  vortices  [22]. 

The  parabolized  stability  equations,  commonly  abbreviated  to  PSE,  were 
conceived  by  Herbert  and  developed  by  Herbert  &  Bertolotti  [20]  to  incorporate 
nonlinearity  and  the  slow  growth  of  the  boundary  layer  into  the  boundary-layer 
stability  computations.  The  results  were  found  to  agree  with  those  of  full  DNS 
simulations  up  to  “spike  stage”,  where  the  complexity  of  the  flow  rapidly  spreads 
beyond  the  resolution  of  the  PSE  code  [11,12,21]. 

We  express  the  disturbance  velocity  field  in  a  series  in  time  (index  /),  in  a: 
(index  n)  and  in  z  (index  fc), 

v{x,y,z,t)  =  £  £  £  Vp(a;,t/,2,<)  (5) 

(=— oo  n=— oo  k=—oo 

where,  p  is  the  wave-vector  (/,  n,  k).  The  velocity  field  of  each  mode  is  parti¬ 
tioned  into 

Vp(x,  y,  Z,  t)  =  Vp(x,  y)  Xp(x,  t).  (6) 

The  function 

Vp(a:,y) 


describes  (in  a  sense  made  specific  below)  the  velocity  profile  of  eax:h  mode  and 
the  function  x  incorporates  the  growth  and  wavelike  part  of  the  mode’s  velocity 


field 


Xp{x,z,t)  =  exp 


s  +  ik^z  —  Hut 


(7) 


where  the  complex  wavenumber  Up  is  composed  of  a  real  growth  rate,  7p,  and 
an  imaginary  wavenumber  ap, 


flp  =  7p  +  iotp 


The  partial  differential  equation  governing  the  velocity  profiles  Vp  is  ob¬ 
tained  by  substituting  the  expansion  (5)  into  equations  (4a)  and  (4b).  For  the 
streamwise  derivatives  we  make  use  of  the  slow  change  with  x  of  the  profiles 
and  growth  rates  with  the  rule 


dx”^ 


flpVp  "f*  fHCtp 


dx 


/  1  \  m— 2  . 

+  if’ 


Xp» 


(8) 


where  ap(x)  =  yp(x:)  ina.  For  m  >  \  the  streamwise  derivatives  of  the 
mean-flow  Vg  are  zero,  in  accordance  with  the  boundary-layer  approximation. 
Similarly,  in  (8)  we  drop  second-  and  higher-order  derivatives  with  respect  to 
X  of  Vp  and  Op.  Performing  harmonic  balance  yields  an  infinite  set  of  coupled 
partial  differential  equations  of  parabolic  type  in  x  of  the  form 

^  E  Q  ,  Vp_r] ,  (9) 

where  the  operators  L,  M,  N  and  Q  contain  derivatives  with  respect  to  y  only. 
The  summation  on  the  r.h.s.  of  (9)  is  is  truncated  to  some  number  (L,  TV,  K)  in 
the  numerical  computations.  Due  to  the  symmetry  in  z  we  only  need  to  solve 
for  modes  with  non-negative  wave  numbers  in  t,  x,  and  z.  Upstream  traveling 
modes  are  not  allowed. 

An  “auxiliary”  condition  is  needed  in  all  PSE  formulations  to  remove  the 
ambiguity  in  (8)  caused  by  the  dependence  of  both  Vp  and  Cp  on  x.  This 
condition  as  also  been  called  a  “normalization”  condition  in  PSE  literature.  In 
our  particular  case  we  employ  the  condition, 

(10) 

where  f  denotes  the  complex  conjugate.  Equations  (9)  and  (10)  form  a  complete 
set  for  the  unknows  Vp  and  Op. 
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We  transfer  the  zero-slip  wall  boundary  condition  to  y  =  0  via  a  Taylor 
series  expansion  about  y  =  0,  and  stop  the  expansion  at  terms  linear  in  \H\ 
since  this  quantity  is  assumed  small,  although  including  higher  orders  in  \H\ 
can  be  done  in  a  straight  forward  fashion.  Performing  harmonic  balance  yields 
the  boundary  conditions  satisfied  by  Vp  for  each  p, 


Vp(x,0) 


-1 

Ap(x) 


avB(x,o) 

**  dy 


(11) 


where 


The  xPSE  transition  analysis  tool-kit  has  been  employed  for  the  computa¬ 
tions.  The  partial  differential  equations  (9)  are  transformed  into  algebraic 
form  by  use  of  a  multi-domain  spectral  collocation  technique  in  y,  and  a  fi¬ 
nite  difference  discretization  in  x.  Five  domains  are  used  in  y,  with  limits 
at  [0,1.6],  [1.6,4],  [4,14],  [14,56],  [56,96]  (recall  the  reference  length  is  6”  at 
R  =  1000),  and  u  and  v  are  approximated  by  17  and  19  Chebychev  polynomials, 
respectively,  in  each  domain.  Asymptotic  boundary  conditions  are  imposed  at 
y  =  96.  For  each  mode  p,  we  construct  the  vector  of  unknows 


Xp  =  {un,u„,a}p 

composed  of  the  Chebychev  coefficients  for  u  and  v,  and  the  complex  wavenum¬ 
ber  a  for  the  mode,  at  the  new  marching  location.  The  nonlinear  algebraic 
system  is  solved  iteratively  by  lagging  the  nonlinear  terms  one  iterate,  and  us¬ 
ing  a  Newton  method  to  solve  the  linear  problem 


VxF  •  Axp  = -F(xp) -b  NLp  (12) 

for  each  mode  p.  Here,  the  term  NLp  represents  the  nonlinear  terms,  and  the 
function  F  the  terms  in  (4a)  and  (4b)  that  are  linear  in  u  and  v. 
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10  Appendix  B.  The  self-similar  vortex 

We  approximate  the  vortex  flow  from  the  leading  edge  to  the  downstream  lo¬ 
cation  at  which  we  begin  the  PSE  simulations,  say  Xo,  with  a  self-similar  form. 
Such  a  form  satisfies  continuity  and  the  boundary  conditions  exactly,  but  pro¬ 
duces  a  small  residual  error  in  the  momentum  equations.  All  terms  of  the 
Navier-Stokes  equations  reduce  to  terms  dependent  on  the  self-similar  variable 
7/  except  some  coefficients,  comminj'  <'rom  z  derivatives,  that  depend  on  x  but 
with  a  variation  sufficiently  slow  to  allow  the  term  to  be  approximated  as  con¬ 
stant.  The  curvature  is  set  to  zero. 

The  case  of  constant  (in  x)  riblet  height  and  the  case  of  height  proportional 
to  y/x,  as  well  as  the  case  of  blowing/suction  of  constant  (in  x)  strength  and 
of  decaying  strength  proportional  to  \j\lx,  yield  quasi  self-similar  solutions. 
Herein  we  will  focus  on  the  constant  amplitude  cases.  The  surface  geometry  in 
the  case  of  constant  riblet  height  is  described  by  the  function, 

//(a:,^)  =  W(o,i)  exp[i^z]  -|-  c.c., 

where  VV(o,i)  is  a  real  quantity  that  controls  the  height  of  the  riblets.  In  the 
case  of  blowing/suction  the  wall  velocity  is  given  by  the  function, 

K(ar,^r)  =  5(0  1)  expli^z]  -f-  c.c.. 


where  «S(o,i)  is  a  real  quantity  that  controls  the  strength  of  the  suction. 
For  the  ribbed  geometry  we  seek  steady  solutions  of  the  form. 


u{x,y,z) 

=  +  c.c. 

(13) 

v{x,y,z) 

=  u,e//(»/)e’^^-  +  c.c. 

X 

(14) 

w{x,y,z) 

=  C.C. 

X  \  X 

(13) 

where  7/  =  y-^xjxQ  is  the  Blasius  similarity  variable.  For  the  suction  case  we 
seek  steady  solutions  of  the  form. 


u{x,y,z)  =  u,e7/(j7)e‘^*-yiii^  -|-  c.c.  (16) 

v{x,y,z)  =  i;,e//(»7)e‘^*  -k  c.c.  (17) 

w{x,y,z)  =  +  c.c.  (18) 
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Reference  length  is  Xq!. 

Equation  (la)  and  equation  (lb)  become  ordinary  differential  equations  with 
independent  variable  rj.  These  equations  are,  respectively, 

aoUaelf  +  0.\^aelS  +  + 

hoVaelf  -f  +  ^>3^/ +  640^  =  0  (19) 

CoU,e//  +  +  C2te“e,/  +  doVatlJ  =  0  (20) 


where  the  primes  denote  differentiation  with  respect  to  r].  The  formulation  (?) 
for  the  constant  riblet  height  yields, 


ao  =  [5/"  +  7,/"'  +  I,*/”  +  C(/  -  -  ,V")  ]/4 

co  =  riri2  +  r/2-C 

ai  =  + vn/2 

C.  =  //2 

a2  =  [v^r-vr+f]/^ 

C2  =  1 

6o  =  -vrn  -  V"'  +  C{C  -f'  +  t?/72) 

do  =  -r 

6,  =  -T,r  -  /72  -  Cfl2 

62  =  2/'  -  f)r/2  -  2C 

63  =  7/2 

64  =  1 

while  the  formulation  (?)  for  the  constant  suction  yields, 


«o  =  (/" + 3vf" + + c(/  -  if  -  i‘n  1/4 

CO  =  7,772  -  772 

01  =  r,{r + 

c.  =  //2 

a2  =  wr-vf'  +  f]l^ 

C2  =  1 

ho  =  -1-  r  +  Crif  +  2C^  ]/2 

do  =  -r 

61  =  -^r  -  /72  -  C//2 

62  =  /'  -  rtf  12  -  2C 

h  =  f/2 

II 

where  /(i/)  is  the  Blasius  stream -function  variable,  and  the  symbol 


xo 


represents  the  coefficient  that  is  dependent  on  x.  The  no-slip  boundary  condi¬ 
tions  on  the  surface  of  the  undulated  wall  are  transferred  to  y  =  0  via  a  Taylor 
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series  expansion.  These  conditions  are. 


«..,/(<>)  =  -H'(o.i)/"(0)  (21) 

»„//(0)  =  0  (22) 

»;,,,(0)  =  0  (23) 

while  the  boundary  conditions  for  the  suction  case  can  be  directly  obtained, 

«.w/(0)  =  0  (24) 

v„l/(0)  =  W(„,,|  (25) 

>';„/(0)  =  0  (26) 


The  condition  on  follows  from  continuity.  The  coefficient  C  prevents  the 
system  of  equations  from  being  truly  self-similar.  As  an  approximation,  we 
neglect  the  x/xq  dependence  in  C  to  arrive  at  a  system  of  ordinary  differential 
equations.  This  approximation  is  acceptable  in  view  of  the  small  values  of  d  of 
interest  -  typically  less  than  0.1  (nondimentionalized  with  6  at  the  PSE  starting 
location  Xg).  Thus,  over  the  region  0  <  x  <  Xo>  the  magnitude  of  C  is  less  than 
<  0.01,  which  is  small  when  compared  to  the  order  0(1)  coefficients  in  (2a) 
and  (2b). 
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